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Abstract 



We built a Low Energy Deuteron Polarimeter (LDP) which mea- 
sures the spin-polarization of deuteron beams in the energy range of 
25 to 80 keV. The LDP works by measuring azimuthal asymmetries in 
the D(d, n) 3 He reaction at 9^ = 90° and comparing them to analyzing 
powers. We built this polarimeter for two reasons. Firstly, to cross- 
check other polarimeters at KVI (LSP and IBP). Secondly to test the 
LDP itself. We were able to use the LDP as a vector polarimeter, but 
not as a tensor polarimeter because of an uncertain tensor analyzing 
power calibration. 




Figure 1: Cartoon depicting polarization. The top panel represents a polarized 
beam, where spins are aligned. The bottom panel represents an unpolarized beam, 
where spins are randomly aligned. 
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1 Introduction 



Most people know that light can be polarized. This means that light oscil- 
lates in a preferred plane. All nuclear physicists know that 'particle beams' 
can be 'spin-polarized'. This means that particles making up a beam have 
their spins aligned in a preferred direction as is illustrated in Fig. 1 (page 
2). 

The linear polarization of light can be determined, quite easily, with a 
polarizing filter. An analogous device, called the Stern-Gerlach filter, can 
determine the polarization of beams of neutral particles by splitting the 
beam according to its spin content. Here the degree of polarization would 
be the excess in quantity of one spin-state over another spin-state. 

Most beams in nuclear physics are made of charged particles such as 
electrons, protons, nuclei, or as in this experiment, deuterons. A deuteron 
is an atomic nucleus consisting of a proton and a neutron. Only beams of 
charged particles can be manipulated (accelerated, steered, focused) by the 
electric and magnetic fields of the machinery in accelerator facilities, such as 
in KVI 1 . Neutral beams are weakly affected by electric and magnetic fields. 
Measuring the polarization of beams of charged particles, as opposed to 
beams of neutral particles, is more difficult because polarization filters such 
as the polaroid or the Stern-Gerlach filters do not exist for charged particles. 
Mott and Pauli (HE], therefore many textbooks in Quantum Mechanics, 
claim that a spin filter for charged particles is theoretically impossible. Yet 
it was recently claimed (Hj that such a device is possible under particular 
conditions. Until such a device is built, or proved unfeasible, the polar- 
ization of charged particle beams is determined by scattering experiments. 
Measuring polarization may seem to be a trivial measurement. But since it 
is a fully- fledged scattering experiment it involves beam, target, detectors, 
and electronics, therefore takes months to accomplish. 

1.1 Nomenclature 

Particle beam polarimetry borrowed nomenclature from optical polarimetry. 
This is not surprising in view of the similarities between spin polarization 
and optical polarization. 

Physics textbooks often describe an experiment with two polarizing fil- 
ters at an angle to each other; the first filter is called the 'polarizer' and the 
second is called 'analyzer'. The wording is such because the first filter can 
polarize (normally unpolarized) light, while the second filter can analyze the 

1 'Kernfysisch Versneller Instituut', The Dutch Nuclear Physics Accelerator Institute 
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strength and direction of the polarization. In nuclear physics polarimetry, 
a reaction of the type 2 A(fe, c)D is called a 'polarization experiment'. The 
vector stands for a polarized specie. In this reaction an unpolarized beam 
b, creates a polarized ejectile c. A reaction of the type A(6, c)D is called 
'analyzing power experiment' because a beam, b, with non-zero polarization 
creates an asymmetry in the ejectile's distribution. The word 'analyzer' of 
optics became 'analyzing power' in nuclear physics. In optics any analyzer 
is as good as the next one, while in nuclear physics different reactions have 
different analyzing efficiencies, hence the word analyzing power. 

1.2 Why Make Polarized Beams? 

Polarized beams or targets are used in nuclear physics to extract observables 
such as analyzing powers and spin-transfer coefficients. These measured ob- 
servables can be compared to predictions of theoretical models to study, for 
example, the three-body force [IJEj, or the spin-terms of the nucleon-nucleon 
potential. 

One application of polarized beams lies in the possibility Ej of using 
polarized deuteron beams to control the reaction rates in fusion reactors and 
to reduce the amount of unwanted neutrons. 

1.3 Our Motivation 

Our motivation to measure polarization stems from the fact that two po- 
larimeters at different beam energies measure different polarizations. One 
polarimeter, called the In-Beam Polarimeter (IBP), measures less polariza- 
tion than another polarimeter, the Lamb-Shift Polarimeter (LSP). By build- 
ing a third polarimeter, this Low-Energy Deuteron Polarimeter (LDP), and 
comparing its results with the other two polarimeters, we wish to determine 
whether one of the polarimeter is inaccurate or whether the polarization 
changes between polarimeters. For beams of polarized protons, the LSP 
routinely measures polarization between 80-90 ± 1% (of the theoretical max- 
imum), while the IBP measures 70-75 ± 4%. Since laboratories, other than 
KVI, are not able to produce polarizations reaching 90%, a healthy skepti- 
cism exists in the LSP's reading of 90%. This new polarimeter, the LDP, 
which measures essentially the same beam as the LSP (in terms of beam en- 
ergy, current and location), can determine whether the LSP systematically 
overestimates the polarization or not. If the LSP and LDP agree on the 
polarization, one can conclude that there are polarization transformations. 

2 The notation A(6, c)D represents a reaction where A is the target, b the beam (projec- 
tile), c and D the observed (ejectile) and unobserved (recoil) products, respectively. This 
notation specifies more than A + b — > c + D can. 
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While if the LDP agrees with the IBP, we conclude that the LSP overesti- 
mates polarization due to some unknown systematic uncertainty. 

Another motivation is that, a deuteron polarimeter based on the D(d, n) 3 He 
reaction has not been published before. Building the LDP is a feasibility 
test for a new polarimeter. The standard reactions, for deuteron polarime- 
try at these energies, are 3 H(d, n) 4 He and D(d,p) 3 H. The advantage of our 
reaction over these two is, respectively, that it does not involve handling of 
tritium, and does need to have detectors inside the scattering chamber. 

1.4 Concept Behind Polarimetry 

By impinging a beam on a target and observing the distribution of scattered 
particles in space, one can deduce the polarization of the beam. 

An unpolarized beam scatters particles with an isotropic (azimuthal) 
distribution, while a polarized beam scatters particles with a non-isotropic 
(azimuthal) distribution. The latter case is illustrated in Fig. [21 




Figure 2: A scattering experiment. A polarized beam impinges on a target. The 
ejectiles scatter preferentially in some direction. Two detectors measure an asym- 
metry in the number of ejectiles. The figure also shows to which parts of the ex- 
periment the quantities P (polarization), A (analyzing power), and e (asymmetry) 
are associated with. 

The degree of non-isotropy in scattering is quantified by the 'asymmetry' 
that detectors measure, and is proportional to the beam polarization with 
a proportionality constant given by the analyzing power. Mathematically, 

e = PA. 
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This equation describes how polarimeters based on scattering, such as 
the LDP and IBP, work. In our experiment P is the quantity which we solve 
for in terms of e and A. e is measured experimentally by the detectors. A 
is a constant that can be calculated or measured by other experiments. 
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2 Theoretical Background 



2.1 Nuclear Potential 

The nuclear potential is not yet fully understood, and is still a subject of 
research and lamentation. An illustrative, but incomplete, form of a two- 
nucleon nuclear potential is 

V = V (r) + V s (r) s[ ■ ? 2 + V SL (r) S • L 

+ Mr) (^(sl • r){sh ■ r) - si • s|) + ^ (1) 

Terms on the right side of the above equation are called central, spin- 
spin, spin-orbit, and tensor. The central term depends only on the distance 
separating the two nucleons. Another example of a central potential is 
gravitational attraction. The spin-spin term, s~i -s~2, stems from the magnetic 
interaction of the spins of the two nucleons. The spin-orbit term, S-L, stems 
from the interaction between the total spin of both nucleons and the angular 
momentum defined by their relative motion. An example of a spin-orbit 
interaction is the 'fine structure' in atomic physics; where the degeneracy for 
states of equal L is lifted. An analogy of a tensor behavior is the interaction 
of two bar magnets. Two bar magnets with parallel orientations placed 
alongside each other (like sardines in a can) will repel each other. Two bar 
magnets with parallel orientation placed along a line (like sardines chasing 
each other) will attract. The tensor term of the nuclear potential exhibits 
the same angle-dependent behavior. 

2.2 Spin-Orbit Term 

The spin-orbit term of the nuclear potential is responsible for the azimuthal 
distributions of scattered particles. Spin-spin and tensor terms contribute 
to other polarization phenomena not addressed is this experiment. In case 
Vsl{t) or S-L is zero in Eq. (J^), one would observe a flat distribution of 
ejectiles in the azimuthal angle. 

Assume, without loss of generality, that Vsl(t) > 0. When S and L 
are parallel, the spin-orbit term is positive, and therefore decreases the at- 
tractive negative nuclear potential. This is illustrated in Fig.|2]a, where the 
trajectory of a particle in a collision gets less deviated in the presence of 
parallel spin and orbit angular momentum. Conversely, if S and L are an- 
tiparallel, then the potential will be more attractive, and the particle will 
scatter more towards the target nucleus as illustrated in Fig.|3]b. In both 
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with spin-orbit 
without spin-orbit 




a) b) 



Figure 3: Scattering of a particle with spin S on a spinless target. The collision 
defines the angular momentum L. In a) S and L parallel, reducing the attractive 
nuclear force. In b) S and L are antiparallel, increasing the attractive nuclear force. 
Without the spin-orbit term trajectories a) and b) would be indistinguishable. 

cases particles are scattered preferentially towards the right direction. 
2.3 Polarization Formalism 

Every deuteron has spin. Furthermore, this spin must be aligned in one of 
three ways with respect to a quantization axes. The spin can be aligned 
parallel (mj = 1), anti-parallel (mi = —1), or perpendicular (mi = 0) to 
this quantization axis. 

The polarization formalism given below is based on articles by Ohlsen 
[SJ |H] and describes the relation between all quantities involved in this 
project. It describes how quantum mechanics, experimental and theoret- 
ical nuclear physics meet. This formalism allows one to derive few simple 
equations (Eq. (JTT|) & Eq. (|12[l) that are used to determine beam polarization 
in terms of known analyzing powers and measured experimental asymme- 
tries. 

Although deuterons are spin-1 particles, I will start the polarization for- 
malism for spin-^ particles then extend it to spin-1 particles. The reason is 
that, the polarization formalism for spin-| particles contains all the essen- 
tial features of the formalism with much less mathematics. Therefore it is a 
good starting point to visualize. 
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X 




r'Z 



(spin quantization axis) 



^■z (beam direction) 



Figure 4: Coordinate system used to describe scattering with a polarized beam. 

2.3.1 Coordinate System 

The coordinate system most often used to describe polarization experiments 
is called the 'Madison Convention' ^U] and is shown in Fig. 2] This coor- 
dinate system incorporates both beam and scattering parameters into one 
coordinate system. The direction of the z-axis is parallel to the momentum 
of the incoming beam, fej n . The y-axis is along fej n x k ou t, where k ou t is the 
direction of the outgoing ejectile. The y-axis, therefore, is perpendicular to 
the scattering plane. The x-axis is left to form a right-handed system with 
the y and z axes. The angle between z and Z is given by 0. The angle <j> is 
between y and the projection of Z onto the x-y plane. Scattering to the the 
left, right, up and down with respect to the quantization axis correspond to 
0= 0, 180, 270 and 90° respectively. 

Note that z (lower-case) is the beam direction, while Z (upper-case) is 
the quantization axis direction. One must distinguish between p z and Pz- 
The former is the component of polarization along the the beam direction, 
while the latter is the degree of polarization along the quantization axis, the 
quantity of interest in this experiment. 

2.3.2 Spin-i Particles 

A spin- \ particle can be represented by a Pauli spinor, 



where a\ is the probability amplitude of finding the particle in the state 



called spin-up, and a>2 is the probability amplitude of finding the 
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particle in the spin-down state f ^ J . 



The spin state of an ensemble of N such particles can be represented by 
a set of Pauli spinors, 

where n runs through all particles. 

If the order of spins in an ensemble is not important, but only the aver- 
age spin, then a beam can be described by a density matrix, 



v^/V !„(»*) |2 
Z^n=l l u l I 
v^Af _(n) _(")* 
l^n=\ a 2 a l 



TV 

n=l 



1 "2 
u( n )|2 



(3) 



A density matrix fully characterizes the polarization (magnitude and di- 
rection) of particle beams. 

Here are some relations, 



ft = Tr{pa x ) = ^Y,2Re{a^at ) *), 

n=l 
1 N 

Py = Tr(p* v ) = ±^2Im(Jp<g>), (4) 

71=1 
1 ^ 

p 2 = Tr(^.) = l^(|4^| 2 -|a^| 2 ). 

n=l 

One can show [Sj that the density matrix can be written as a linear 
combination of Pauli spin operators, 



P 2 



9 (J + (5) 

where: 

• pj 's are the components of polarization in the x, y, and z directions; 

• / is the unit, mat rix ( J 
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<jj 's are the Pauli spin operators; 





' 






cr x = ( 




i 


; CTy = f 



o -i\ MO 
• o ; ° z = [ -1 



A nuclear reaction can transform the spin state of particles. Therefore, 
the spinor of an outgoing particle X/ is related to the spinor of an incoming 
particle Xi by a transformation, M, 

Xf = M X i- 

M is a 2 x 2 matrix whose elements are functions of energy and angle. 

The density matrix describing the incoming beam pi can be written in 
terms of spinors as 

N 

« = Exi B) [xS B) i t , 

n=l 

and for the outgoing beam, 

N 

>)k>)lt 



n , — V J n > N/ (n 



n=l 

The density matrix is transformed by a reaction as 



Pf = MpiMl (6) 



The cross section of the reaction can be given by, 

^.fl-^-^ftf). (7) 



If the beam is unpolarized 

P i oc 



1 
1 



and if the density matrix is normalized to unity 

1/10 



Pi ~ 2 V 1 
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then Eq. (JJJ) reduces to 

do = —Tr(M Aft) (unpolarized beam). 

Applying Eq. © to Eq. © one gets, 

1 1 3 

p f = -M Aft + - Af^Aft. 

i=i 

Taking the trace yields 

3 

a(9, <P) = Tr{ Pf ) = a {l + M e ))> ( 8 ) 

i=i 

where 

_ Tr(MajM^) 
j ~ Tr(AfAft) J-^y^- 

Af is the same 2x2 matrix that was used to transform the spin state in 
Eq. ((2) and density matrix in Eq. (jfjj). A/'s are the analyzing powers of the 
reaction. The analyzing powers, like cross section, are a property of nuclear 
reactions. They can be calculated from theory or measured experimentally. 

Parity and time reversal arguments [S] reduce Eq. (jHJ) for two-body reac- 
tions to 

a(9,^=a (l+p y A y (d)). (9) 

This implies that reactions are only sensitive to polarization along the y- 
axis. Polarization along the x, or z direction do not affect the cross section. 

Since p y = Pz • y = Pz cos 4>, Eq. Q can be written as 

(7(0,0) = a (l + P z cos 4>A y (9)), (10) 

which is plotted in Fig. [21. 

Scattering to the left of the quantization axis {(j) = 0) has a cross section 

a L = ao(l + P z A y ). 
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Figure 5: Azimuthal angle dependence of the cross section for Pz = 1 and A y = 
0.5 . In this example there are more particles going to the left direction than to the 
right direction. Detectors placed in these directions measure an asymmetry. 



A detector placed in this direction will measure a count proportional to this 
cross section 

L oc ol- 

Scattering to the right of the quantization axis (</> = 180) would have a 
cross section of, 

or = o"o(l - PzA y ). 

A detector placed in this direction will measure a count proportional to 

R oc Or. 

Define the asymmetry, e±, of the reaction as, 

L-R 

£i = lTr- 

The asymmetry is, therefore, the difference over the sum of the number of 
particles detected in the left and right detectors. 

Scattering in a direction parallel to the quantization axis has the same 
cross section as scattering in an anti-parallel direction. Therefore two detec- 
tors placed above and below the scattering will measure the same number 
of ejectiles, resulting in no asymmetry. 
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It can be seen that 

P Z Ay = e 1 . (11) 

This equation relates the polarization, analyzing power, and asymmetry for 
a spin-i beam. We introduced this equation in Sec. 11.41 and now present it 
with subscripts reflecting some of the geometry behind the scattering. 

2.3.3 Spin-1 Particles 

The cross section of a reaction using a spin-1 polarized beam is 

/ 2 3 
a{0,<j>) = (7 (1 + ^PyA y {9) + -p xz A xz (6) 

+ -£PxxAxx{6) ~ \pyyA yy (0) + \pz Z A zz (6)*j . 

This equation is the spin-1 analogy of Eq. Q. Yet it contains noticeably 
more terms! 

By expressing the polarization in terms of the coordinates (Pz, Pzz, P, 4>) 
instead of (p y ,p xz ,p xx ,p yy ,p zz ) we get, 

a(6, 4>) = do ^1 + —PzAy(0) cos sin p — PzzA xz (6) sin j3 cos /3sin</> 
- \Pzz(A xx (9) - A yy (6)) sin 2 (3 cos 2<j) + ^P zz A zz (e)(3cos 2 (3 - l)). 

In the spin-^ case, only one polarization (Pz) and one analyzing power 
(Ay) enter the equations, while for spin-1 beams, two polarizations (Pz, 
Pzz) and four analyzing powers (A y , A xz , A xx — Ay y , A zz ) enter the equa- 
tions. This complexity arises because a spin-1 beam has three spin substates, 
while a spin-^ beam only two (spin up and spin down). 

It is worth mentioning that quantities with one index are called vectors, 
while quantities with two indices are called tensors. For example, Pz is 
called 'vector polarization', and A zz is one of the three 'tensor analyzing 
powers'. 

Scattering to the left, right, up, and down (4> = 0°, 180°, 270°, 90°) 
directions with respect to the quantization axis have cross sections, 
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°L = <7 (1 + ^PzA y (9) sin/? + -P zz (A yy (0) sin 2 + A zz cos /?)), 



= a (l- -Pz A y (0) sin 13 +-Pzz(A yy (0) sin 2 /3 + A zz cos 2 /?)), 



ac/ = CT O (l + PzzA ;EZ (0)sm^cos^ + -Pzz(^(^) sin 2 /? + ,4 Z2 cos 2 /?)), 
ct d = a (l + PzzA xz {9) sin (3 + -P Z z(A yy {9) sin 2 /? cos/? + A Z2 cos 2 /?)). 



Detectors placed in these scattering directions will measure a count (L, P, 
{/, D) proportional to the cross sections, 



L oc (Jl 

R OC £TR 
[/ OC £7(7 
-D OC (T£). 

One defines the five asymmetries ei, €2, £3, £4, €5, by the following equations 



e 2 



L- R _ %P z sinf3A y 



L + R \ + \P zz [sin 2 /3A W + cos 2 /3A ZZ ] 

U — D = P zz sin (3 cos(3A xz 

U + D ~ 1 + ip zz [sin 2 f3A xx + cos 2 /JA 22 ] 



2(L - g) fPzsin^Ay 



£4 



£5 



L + P + + D 1 + ±P zz [3(cos 2 /3-l)A 22 ] 

2(t/ - g) = P zz sin/3 cos ft 
L + R + U + D l+ ip zz [ 3 (cos 2 /?-l)A 2Z ] 

(L + P) - (17 + £>) -\Pzz sin 2 /? (A^ - A yy ) 



L + R + U + D 1 + ip zz [3(cos 2 /? - 1)A Z 



which relate the asymmetries, analyzing powers, and polarizations for spin-1 
beams. These are the generalization from spin-| (Eq. IjlljO to spin-1 parti- 
cles. 



2.4 d + d Reactions 

What happens when a deuteron beam strikes a deuteron target? Many 
things happen, so let's restrict the question to: What nuclear reactions are 
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Reaction type 


Reaction Q-value (MeV) 


Elastic scattering 1 


B(d, d)D 





One-particle rearrangement 


D(d, n) 3 He 


3.27 




D(d,p) 3 R 


4.03 


Radiative capture 


D(d, 7 ) 4 He 


23.8 


Breakup reaction 2 


B(d, n)pB 


-2.2 



Table 1: d + d nuclear reactions. 



possible when a deuteron beam strikes a deuteron target? Tabled lists all 
known reactions. 

A number of reaction listed in Table ^ take place in the center of the 
sun, and possibly in future fusion reactors. Coincident ally, solar and reactor 
plasmas are in the same energy range as in this experiment. Yet this setup 
is used to measure polarization, and not to study solar plasmas. 

The cross section of the B(d, n) 3 He reaction highly depends on the inci- 
dent deuteron energy, as is shown in Fig. El A simple quantum mechanical 
effect can model this dependence satisfactorily, as is explained below. 

The potential between two nuclei is attractive at short distances (a few 
fm) because of the strong force, and repulsive at large distances (greater than 
a few fm) because of the coulomb force. A deuteron approaching another 
deuteron sees a potential barrier with a height of about 1 MeV. How can a 
deuteron with kinetic energy in keV range get through the barrier to fuse 
with the other deuteron? The answer is quantum mechanical tunnelling! 
The tunnelling probability, hence the cross section, is proportional to 

a{E) oc e^l (13) 

This dependence is called the 'Gamov factor'. Fitting it to published data 
sets [TU of D(d,n) 3 He cross sections yields an empirical relation, 

a{E)=ae^E (14) 

with a = 22 ± 4 mb and b = 30 ± 2 VkeV. 

In elastic scattering total kinetic energy is conserved. For example, billiard ball col- 
lisions. All particles in D(d, d)D are deuterons. In this notation capital letters represent 
an atom or molecule, while lower-case letters represent a nucleus. D is called deuterium, 
while d is called deuteron. These particles differ only by one electron. 

2 This is the only reaction that does not take place at our beam energy. Beam energy 
needs to exceed the Q-value for the reaction to occur, due to conservation of energy. 
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2.5 Stopping Power 

The basis of any understanding in experimental nuclear and particle physics 
depends on the understanding of the 'passage of radiation though matter'. 
Instances of 'passage of radiation through matter' are; an alpha beam im- 
pinging on a gold foil (Rutherford's famous experiment), a neutron deposit- 
ing energy in a detector (as in this experiment), and ion radiotherapy (where 
ion beams destroy tumors). 

As charged particles pass through matter they lose energy and are de- 
flected. This is not surprising in light of the many imaginable ways in which 
particle and matter can interact. The most important phenomena that con- 
tribute to the net process are 

1. Inelastic collisions with atomic electrons 

2. Elastic scattering from nuclei 

3. Inelastic nuclear reactions 

4. Cherenkov radiation 

5. Bremsstrahlung 

6. etc . . . 

A formula which models these phenomena is the well-known Bethe-Bloch 
formula ^21- This semi-empirical formula gives the energy loss (§§) of 
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various charged particles though various materials. With the energy loss — 
also called stopping power — one can calculate the range of particles and 
total energy deposited in matter. These quantities are vital for detector 
and safety consideration, since you want to know where your particles are 
going and with how much energy! Although the Bethe-Bloch formula for 
stopping power is extensively used in nuclear physics, it is only valid for 
particles with energies greater than 1 MeV/nucleon. At the low energies of 
this experiment another model and formula for stopping power is given by 
Lindhard 13 , which states that the stopping power is proportional to the 
beam velocity 

dE i— 

— — ex v oc v E. (15) 

ax 

Anderson Jl] has compiled empirical formulas for stopping powers of 
various beam and target combinations, by fitting data with the Bethe-Bloch 
and Lindhard formulas. The stopping power of a deuteron beam on a C2D4 
target is given 3 by 

dE 

— = 14.48 E 0A5 ([E] = keV and [x] = /im). (16) 

ax 

Calculation of the mean penetration depth of the beam into the target 
yields 0.4 to 0.7 /im for Ed = 25 to 80 keV. This is the distance within which 
the mean beam energy decreases to zero. Since many relevant parameters — 
such as cross section and analyzing power — are energy dependent, their av- 
erage value has to be calculated. The cross section is energy dependent, and 
since there is an energy loss of the beam through the target, therefore the 
average cross section must be calculated. 

The Bethe-Bloch Formula does not apply to neutral particles such as 
photons and neutrons. Neutral particles have a larger range through matter 
than charged particles. Particles created from a d + d reaction, at these 
energies, include neutrons (n), protons (p), tritons ( 3 H), helium-3 ( 3 He), 
helium ( 4 He), and 7-rays. 7-rays are also observed from the de-excitation 
of nuclei after having absorbed neutrons. Therefore, wherever one observes 
neutrons, one is likely to observe 7-rays as well. The charged particles (p, 
3 H, 3 He, and 4 He) do not make it out of the vacuum chamber. They are 
stopped in the glass beam tube, metal target holder, or the target itself due 
to their small range of order order /im to mm. This tiny range is due to 
the high energy loss of low energy charged particles passing through matter. 
Neutral particles, on the other hand, such as neutrons and 7-rays, do not 
interact as much with material and are able to exit the beam tube and reach 

3 Use Bragg's Rule and the the assumption that the energy loss of deuterons through 
deuterium is equal to that of protons through hydrogen. 
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the detectors or go beyond them. This is why the D(d, n) 3 He reaction was 
chosen, because it could be isolated from other reactions also occurring. 
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3 Experimental Set-Up 

The different devices of our experiment are sketched in Fig.[7| 



POLIS 




LSP 




LDP 




AGOR 




IBP 











Ion Source Polarimeter Polarimeter Accelerator Polarimeter 



Figure 7: Schematics of the beam line we used. Polarized beams are created by 
POLIS, the beam polarization is measured with the LSP and the LDP, the beam 
is accelerated, and finaly the beam polarization is measured with the IBP. 

We briefly describe POLIS. We describe at length the components of the 
LDP. The LSP gMUM, AGOR [J7|, and the IBP [TBI are described elsewhere. 

3.1 Polarized Ion Source 

Our Polarized Ion Source ^Hj (POLIS) provides beams of polarized protons 
or deuterons. Proton beams can be vector polarized, while deuteron beams 
can be vector and/or tensor polarized. 

Before the invention of polarized ion sources, such as POLIS, polarized 
beams were produced by using the scattered ejectiles of a reaction. These 
ejectiles were partly polarized, and were used themselves as a beam for an- 
other experiment. These 'double scattering' experiments were plagued by 
low polarization and intensity. 

The production of a polarized deuteron beam from POLIS is outlined as 
follows, 

„ dissociation „ transitions ionization — f 

D 2 > D ► D* ► d 

First, deuterium molecules from a gas cylinder are dissociated into deu- 
terium atoms and collimated into a beam. Then, electromagnetic transitions 
polarize the atom by populating some of its hyperfme states. The polarized 
atoms are ionized leaving only a beam of polarized deuteron nuclei. This 
beam is then accelerated to desired energies, and steered to a particular 
experiment. 

Each hyperfme state of an atom represents a particular alignment of the 
spin of the nucleus and electrons. The electromagnetic transitions between 
hyperfme states are produced by applying RF fields to atoms present in a 
magnetic field, as depicted in Fig.EJ POLIS has four transitions units called 
Weak Field, Medium Field, Strong Field I, and Strong Field II. Different 
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Figure 8: Breit-Rabi diagram of a deuterium atom in the ground state. The spin 
substate of the nucleus (d) and its electron (e) are listed on the right of the figure. 





POLIS 


Max. 


Theoretical 


Beam 




Transition 


Polarization 


Polarization 




Units 


Pz 


Pzz 






WF 


-2/3 







Positive Vector 




SF. I + SF. II 


2/3 







Negative Vector 


MF 


+ sextupole + SF. I 







1 


Positive Tensor 


MF - 


f sextupole + SF. II 







-2 


Negative Tensor 



Table 2: POLIS transition units with the beam polarization they generate. Other 
combination of transition units are possible yielding beams with both non-zero Pz 
and P Z z- 



combination of these transition units populates different hyperflne states. 
Populating a particular hyperflne states is selecting a spin substate of the 
nucleus and electrons. The spin substate of electrons is irrelevant since elec- 
trons are stripped from the atom. The spin substate of the nuclei constitute 
the beam polarization. 

3.2 Low Energy Deuteron Polarimeter 

The LDP consists of a target and four detectors as shown in Fig.|§J A polar- 
ized deuteron beam from POLIS impinges on a target containing deuteron 
in the form of C2D4 to reacts as D(d, n) 3 He sending neutrons in space. 

3.2.1 Target 

A deuterated-polyethylene (C2D4) film was used as a deuterium target. It 
has the same chemical properties as polyethylene (i.e. the hydrocarbon 
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Figure 9: Photograph of the LDP showing beam tube, target holder and four 
neutron detectors. Also shown are the directions of the beam, the quantization 
axis, and left, right, up, and down scattering. 

C2H4), except that it contains deuterium nuclei (deuterons) instead of hy- 
drogen nuclei (protons). Although carbon is present in the target material, 
it does not produce neutrons in a nuclear reaction because the Q-value of 
the 12 C(d, ra) 13 N reaction (Q = —281 keV) is higher than our beam en- 
ergy. Therefore neutrons detected originate exclusively from the D(d, n) 3 He 
reaction. Alternative targets to C2D4 exist, such as, deuterated-titanium, 
or deuterium gas targets. These were not used because C2D4 targets were 
readily available at the KVI and are more convenient to produce. Target 
thickness was in the order of a few hundred /Ltg/cm 2 , which corresponds to 
a target depth of a few //m. 

At early stages of the experiment, targets consisted of a C2D4 thin film 
held by a rectangular frame. Eventually targets evolved into a C2D4 thin 
film on a round metal backing. There were two problems with the initial 
target design. First, the target would melt under beam heating. The reason 
was that Polyethylene, being a plastic, conducts poorly the energy deposited 
by the beam. The solution was to couple the thin-film to a metallic backing 
acting as a heat sink. Secondly, the sharp edges of the rectangular frame help 
produce unwanted electrical discharge when the target was at High Voltage. 
A round metal holder increased the breakdown voltage. The breakdown 
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voltage 4 was found to be proportional to the pressure in the evacuated beam 
line, the deuteron beam current, and surface conditions on the target. 

3.2.2 Detectors 

The detectors we used were 'liquid organic scintillators' of type NE213. 
These detectors are frequently used for the detection of neutrons. The sig- 
nal produced by these detectors depends on the type of particle entering 
them. 



3.2.3 Pulse Shape Discrimination 

As determined earlier (Sec. l2,5l and 12,4) ). only neutrons and 7-rays will reach 
the detectors. One wants to distinguish between detected neutrons and 7- 
rays because one needs to measure the asymmetry originating from a single 
reaction, and not two reactions. In our case neutrons should be counted, 
while 7-rays should be rejected. A method called Pulse Shape Discrimina- 
tion |12l I2U| allows different particles to be distinguished based on the signal 
they produce in detectors. Different particles have different energy loss 
mechanisms inside matter, and so produce sightly different signal shapes. 
Neutrons will deposit their energy more slowly than gamma-rays, therefore 
the signal they create decays more slowly. The signals are illustrated in 
Fig.Eni 

By monitoring the signal shape which the detected particles produce, 
one can identify particles. We monitored the signal shape by integrating two 
copies of the signal with two different time intervals, and taking their ratio. 
The signal and time intervals are illustrated in Fig. llll The ratio of these 
two integrated quantities is proportional to the decay time of the signal, and 
is the basis for particle identification. By plotting the occurrence of signals 
as a function of this ratio we get a pulse shape spectrum as shown in Fig. 1121 

We can also plot the integral of the signal with the long gate versus that 
with the short gate to get a 2-D scatter plot of pulse shapes, as shown in 
Fig.Hl 

3.2.4 Post Acceleration 

Placing the target at a potential is an experimental trick to reduce the mea- 
suring time of our experiment. The cross section of the D(<i, n) 3 He reaction 

4 A high voltage electrode in pressures of 10 -2 to 10~ 3 mbar generates all kinds of 
plasma effects that are beautiful to watch, such as, striations and micro-discharges. 
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TIME [ns] 

Figure 10: Pulse shape differences for neutron and 7-rays. Detector signals are 
shown on the right panel and the integral of the signals on the left panel. 




short gate 



30 ns 
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Figure 11: Detector signal and integration gates as used in our experiment. 
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Long/Short Pulse Component 

Figure 12: Pulse shape spectrum for one detector. Neutrons are well separated 
from a large background of 7-rays. The number of neutrons detected is the number 
of events under the neutron peak. 
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Figure 13: 2-D pulse shape spectrum. Detected neutrons can clearly be separated 
from detected 7-rays. 

was found to depend on the reaction energy (Fig.EJ). By increasing the en- 
ergy of the beam, the cross section, therefore reaction rate, increases. High 
counting rates have the obvious advantage of lower statistical uncertainty 
for a given measuring time. The maximum beam energy that POLIS can 
produce is 30 keV. Placing the target at negative potentials will accelerate 
the positively charged deuteron beam. If the beam is initially at 30 keV 
and the target at —50 kV, then the beam energy as it reaches the target is 
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80 keV. The reaction rate at 80 keV is about an order of magnitude higher 
than at 30 keV! 

Post acceleration equipment is not shown in the photograph of the LDP 
(Fig.lHJ). It consists of nothing more than a high voltage cable connecting 
the target to a high voltage power supply surrounded by safety features 
(insulation and grounding). 

3.2.5 Data Acquisition 

Data acquisition consists of electronic modules, a CAMAC crate, and a 
PC. Electronic modules can manipulate and perform basic operations on 
electronic signals. The detector signal was first split into multiple copies 
by a 'fan-in/fan-out' unit. Two copies of the signal were each put into a 
'constant fraction discriminator' (CFD) unit to generate logical gates. One 
long and one short gate acts as integration windows for the detector signal 
as shown in Fig. 1111 Another copy of the signal was cable delayed before 
being put to a 'charge to digital converter' (QDC) units along with the short 
and long gate. The QDC integrates the voltage of a signal and converts it to 
into a binary string which can be read by a PC. The program DAX, based 
on the CERN Program Library package |21] . was used to write raw data 
into event files (*.ntuple). PAW |22| was used to perform a simple analysis 
of the data such as plotting and counting. 
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4 Data Analysis 



4.1 Low-Energy Deuteron Polarimeter 
4.1.1 Analyzing Powers 

Analyzing powers are a manifestation of the dependence of the reaction 
cross section on spin. The analyzing powers relevant to this experiment are 
almost entirely available in the literature. Becker |23| has measured all an- 
alyzing powers (A y , A xz , A zz , A xx — A yy ) for the D(d, n) 3 He and D(d,p) 3 H 
reactions at the reaction energy of E = 28 keV. Fletcher |24l I25j has mea- 
sured two tensor analyzing powers (A zz , A xx — A yy ) of both reactions at the 
beam energies of Ed = 25, 40, 60, and 80 keV. Tagashi j^H] has measured 
all analyzing powers of the D(d,p) 3 H reaction at Ed = 30, 50, 70, and 90 
keV. Our polarimeter is limited to the energy range of Ed = 25 to 80 keV, 
because that is the energy range at which analyzing powers are presently 
known, also because those are the beam energies available to us. 

A xz and the energy dependence of A y are not published. A y is claimed 
|23| to be energy independent in our energy range. Furthermore, A y , for 
the D(d, p) 3 H reaction, has negligible energy variation as can be seen from 
|26| . An unknown A xz is not a problem for us since we do not use it in our 
analysis. When (3 = 90° (as in our setting of the beam) only ei, £3, and €5 
enter the analysis (Ea.ll2|), 

Analyzing powers are usually reported in the literature as data points 
with fitted curve (Fig.ll4|) as in |23[ I25j . or just fitted curves (Fig. ll5j) as in 
|24l I26j . The fitting functions are Legendre polynomials. 
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Figure 14: Data and fit of A y at 28 keV from Becker [23]. 
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Figure 15: Fit of A zz at various beam energies from Fletcher |25| . 



4.1.2 Effective Analyzing Powers 

In our analysis one cannot directly use the analyzing powers as found in the 
literature. The asymmetries measured by the detectors are the convolution 
of both non-zero analyzing powers and experimental effects, such as: 

1. Finite acceptance of the detectors. Detectors cover a non-zero solid angle. 
Therefore, the analyzing power must be averaged over this solid angle. We 
chose to average A(0) over A8, and disregard A(j> effects because they are 
small at Za6 = 90°. 

2. Change in analyzing power due to the energy loss of the beam through the 
target, A(E). Fig.lTKl shows that A zz is energy dependent. 

3. Change in cross section due to energy loss of the beam though the target, 
a(E). 

4. Energy profile of the beam through the target; E(x). 



These experimental effects give rise to an effective analyzing power given 

by 

A eff = lAeL x A (E^),e) *(E(x))d9dx 
A6 J Ax a[E(x))dx 

The opening angle of each of our detectors is A# = 26°. Ax is the range 
of the beam. E(x) and Ax can be calculated from the stopping power of a 
deuteron beam on a C2D4 target (Eq. ([160 ■ E(0) = E& and E(range) = 0). 
A(E, 9) represents the energy and polar angle dependence of any of the four 
analyzing powers (A y , A xz , A zz , A xx — A yy ). A{6) is given in the literature. 
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E d = 50 keV 


E d = 80 keV 


A°tf 

y 


0.186 ±0.01 


0.186 ±0.01 


zz 


0.31 ±0.03 


0.25 ±0.03 


(A -A ) e M 

\-<*-xx -^yy ) 


-0.98 ±0.1 


-0.83 ±0.1 



Table 3: Effective analyzing power of the LDP for 50 and 80 keV beams. 

A(E) was estimated by fitting a quadratic polynomial through the four en- 
ergies at which analyzing powers are known. 

The effective analyzing powers at two different beam energies are given 
in TableEl 



4.1.3 Instrumental Asymmetries 

Fig. 1161 shows the neutron peaks in each detector coming from the reaction 
with a polarized and unpolarized beam. 

One can observe that the neutron peaks with an unpolarized beam are 
not of the same height, indicating asymmetries! For example, the peak 
of the Right and Up detector approach 2000 neutrons per bin, while the 
Left and Down peaks approach 1000 neutrons. This is due to instrumental 
asymmetries. One should distinguish between reaction asymmetries and in- 
strumental asymmetries. Reaction asymmetries are fundamental and due to 
a non-zero analyzing power of the nuclear reaction. Instrumental asymme- 
tries are due to mismatches in detector response, such as differences between 
detectors in, gain, solid angle, detection efficiency, discriminator threshold 
level, or signal shape. The fact that instrumental asymmetries are equal to 
a factor of 2 here, indicates that the detector responses were not matched 
properly. This is not a problem because one uses unpolarized beams to 
quantify instrumental asymmetries. Unpolarized beams theoretically do not 
create reaction asymmetries. Therefore, any asymmetry measured with un- 
polarized beams is taken to be instrumental asymmetries. Normalizing the 
detector counts obtained with polarized beams, by counts obtained from 
unpolarized beams, leaves only the reaction asymmetries needed to deter- 
mine beam polarization. Therefore, for the right detector, the number of 
neutrons creating reaction asymmetries is, 

„ Rpolarized 

K = ~R ' 

^unpolarized 

And similarly for the other three detectors L, U, and D. 

Until this point, we have covered the effective analyzing powers and the 
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Figure 16: Neutron peak in each of the LDP's detectors. The dashed line stems 
from an unpolarized beam, whereas the solid line stems from a positive vector 
polarized beam. 



normalized asymmetries. Using Eq. (jl2() . one solves a system of equations 
to get the unknowns (Pz, Pzz)- 

4.1.4 Sample Calculation 

As an example of the procedure to calculate the polarization from neutron 
counts, we analyze the data presented in Fig.lTol 



R = 62566/51124 = 1.22380, 

L = 25815/26985 = 0.95664, 

U = 56244/51337 = 1.0955, 

D = 26600/23455 = 1.1340. 

With the asymmetries, 
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ei 

£4 
£5 



L-R 
L + R 
U-D 



= 0.1225 = 



0.0172 



U + D 

2(L - R) 
L+R+U+D 

2{U-D) 
L+R+U+D 
(L + R)- (U + D) 
L+R+U+D 



0.279 P z sin 

1 + \Pzz [-sin 2 /3 0.647 + cos 2 (5 0.312] 
_ 0.0109 P ZZ sin/3 cos (3 

" 1 + \Pzz [sin 2 /? 0.958 + cos 2 p 0.312] 
0.279 P z sin (3 



= 0.1211 



= -0.0174 



-0.0111 



1 + 0.078 P zz [3(cos 2 /3- 1)0.312] 

0.0109 P ZZ sin (3 cos /? 
" 1 + 0.078 P zz [3(cos 2 /?- 1)0.312] 
0.245 P zz sin 2 ^ 



1 + 0.078 P zz [3(cos 2 /3- 1)0.312] 



Fixing /? = 90° and solving for Pzz in £5 yields P^z = —0.05. Solving 
for Pz in €3 or ei using Pzz yields Pz = 0.4. Applying propagation of errors 
to these equations with the statistical uncertainty originating from counting 
uncertainty (AN = y/~N), and the systematic uncertainty originating from 
the uncertainty in analyzing powers given in Tabled yields, AP|z* = 0.001, 

A pstat _ q QQ 3) Ap sys _ Q 005; Ap sys _ Q Q2 Thig measurement Q f 

polarization is plotted in Fig.HTl 
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Figure 17: Measured polarization and maximum theoretical polarization for a 
positive vector polarized beam. The vertical axis is the vector component of po- 
larization, while the horizontal axis is the tensor component. In this plot vector 
polarized beams lie on, or close to, the vertical axis, while tensor polarized beams 
lie on, or close to, the horizontal axis. 
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4.1.5 Low Energy Deuteron Polarimeter 

The Low Energy Deuteron Polarimeter (LDP) data are tabulated in TableE] 
and plotted as squares in Fig.lTHl 



POLIS 




Detector 




Measured 


State 




Count 




Polarization 




R 


L 


U 


D 


stat sys 

Pz ± ± 


stat sys 

Pzz ± ± 


WF 


30703 


43889 


22943 


51843 


0.49 ± 0.0008 ± 0.03 


-0.095 ±0.003 ±0.005 


ST. I + ST. II 


25815 


62566 


26600 


56244 


-0.45 ±0.001 ±0.03 


-0.043 ± 0.003 ± 0.0006 


MF + ST. I 


29379 


53164 


20200 


41951 


-0.030 ±0.001 ±0.005 


0.414 ±0.003 ±0.05 


MF + ST. II 


18011 


33064 


24853 


54641 


-0.008 ± 0.0009 ± 0.002 


-1.1 ±0.004 ±0.1 


Unpolarized 


54312 


100427 


444617 


101696 









Table 4: LDP counts and polarizations. 
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Figure 18: LDP measurements of polarization of four polarized beams. Systematic 
error bars are just visible, while statistical error bars are smaller than the symbol 
size. 



4.2 In-Beam Polarimeter 

The In-Beam Polarimeter (IBP) measures similarly to the LDP. Namely, 
by measuring asymmetries and exploiting known analyzing powers. The 
reaction used to measure asymmetries was H(<i, d)p at = 80 MeV. The 
analyzing powers for this reaction, at this energy, are not reported in the 
literature, so we used calculated [^7] analyzing powers shown in Table |SJ 
The uncertainties in analyzing power (~ 3%) were taken as the variation in 
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analyzing power between various potentials. 

The IBP data are tabulated in Tableland plotted as stars in Fig. 1191 





E d = 80 MeV 
cm = 100° 


Ay 


-0.291 ±0.005 




-0.74 ±0.02 


A A 


-0.107 ±0.003 



Table 5: IBP analyzing powers. 
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Detector 
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State 




Count 




Polarization 




R 


L 


U 


D 


stat sys 

Pz ± ± 


stat sys 

Pzz ± ± 


WF 


29589 


17220 


25004 


24368 


0.478 ± 0.0008 ± 0.01 


-0.443 ±0.006 ±0.001 


ST. I + ST. II 


13180 


20140 


18465 


18085 


-0.643 ± 0.001 ± 0.008 


-0.112 ±0.005 ±0.004 


MF + ST. I 


19407 


17626 


18181 


17932 


-0.131 ±0.001 ±0.003 


0.478 ± 0.007 ± 0.004 


MF + ST. II 


10834 


10257 


13681 


13391 


-0.112 ±0.001 ±0.002 


-1.44 ±0.009 ±0.01 


Unpolarized 


71690 


73594 


75327 


75093 









Table 6: IBP counts and polarizations. 
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Figure 19: IBP measurement of polarization of four polarized beams. Error bars 
are not visible. 
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4.3 Lamb-Shift Polarimeter 



The Lamb-Shift Polarimeter (LSP) uses a different technique, than the LDP 
and IBP, to measure polarization. It measures directly the spin substate dis- 
tribution of a beam, from which the polarization can easily be obtained. 




Figure 20: Spin state distribution of a negative vector polarized beam measured 
by the LSP. 

Let iV_|_ be the population of spins in the spin substate mj = +1. No 
the population of spins in the mi = state, and iV_ for mi = —1. The 
polarization in terms of populations is, 

N + - iV_ 
N+ + N + Nl 
N 

1 ~ 3 N + +N + N_' 

Populations are calculated by fitting three gaussian peaks and a base- 
line to the LSP's spin substate distribution, as is illustrated in Fig. 1201 The 
area under the peaks are equal to the populations N + , No, and N_. The 
uncertainty in fitting parameters (obtained from the fitting program) were 
propagated into uncertainties in polarization. The data shown in Fig. 1201 
yield P zz = -0.121 ±0.004 and P z = -0.51 ±0.01. These uncertainties are 
purely statistical. 

The LSP data are tabulated in Tableland plotted as triangles in Fig. 1211 



Pz = 
Pzz = 
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Pz 


Pzz 


VV -T 
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vj.tju zn u .u± 


u.±U(j zn u.uuo 


ST. I + ST 


II 


0.52 ± 0.01 


-0.141 ± 0.005 


ST. I + ST 


II 


0.510 ±0.01 


-0.144 ±0.004 


ST. I + ST 


II 


0.507 ±0.003 


-0.133 ±0.004 


MF + ST. 


I 


-0.009 ±0.009 


0.68 ±0.07 


MF + ST. 


I 


0±0.07 


0.67 ±0.07 


MF + ST. 


II 


0±0.01 


-1.50 ±0.05 


MF + ST. 


II 


0±0.02 


-1.52 ±0.05 



Table 7: LSP data, taken within a few minutes of each other. 
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Figure 21: LSP measurement of polarization of four polarized beams. 
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5 Results 



The measurements taken by the three polarimeters in Fig. 1181 Fig. 1191 and 
Fig. are superimposed in Fig. 1221 
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Figure 22: Polarization measurements of all three polarimeters. 

These data were collected on the night of 15-16 October 2004. LSP mea- 
surements were taken between 10 and 11 pm, IBP measurements were taken 
between 3 and 5 am, and LDP measurements were taken between 6 and 11 
am. The LSP and LDP disagree on the tensor component of polarization, 
while the LSP and IBP disagree on the vector and tensor components of 
polarization. It is difficult to draw definitive conclusions from this plot for 
three reasons: 

1. The LDP does not fully agree either with the IBP or LSP, not allowing 
one to prove whether the polarization changes between LSP and IBP, 
or if the LSP or IBP have unknown systematic errors. 

2. We do not have a big data set with multiple measurements at this 
moment. 

3. We cannot exclude changes in beam polarization while the polarime- 
ters were being switched. Had we monitored the polarization with 
one polarimeter (LSP) before and after the measurements of other 
polarimeters, we would know whether the beam polarization was con- 
stant over the measurement period. 
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Repeating the experiment would have been the natural way to proceed, 
this would have addressed points 2 and 3. Yet because of time and facility 
constraints, we have to evaluate the present data. 

The IBP data are not reliable for two reasons. Firstly, we do not know 
how successful the theory, which calculates IBP's analyzing powers, is. Sec- 
ondly, the asymmetries we measured depended on how we analyzed the data 
(where we placed cuts). Both of these issues can be addressed with more 
time. Until then we cannot be confident about the present IBP polariza- 
tion measurements and uncertainty in polarizations. Therefore we cannot 
presently answer the question whether there are polarization changes in the 
accelerator. 

One can observe that the polarization measurements of the LSP and 
LDP agree quite well for the Pz component of polarization but not for Pzz- 
One can also see that the LDP measures always less Pzz than the LSP (all 
LDP points are closer to the vertical axis in Fig.l22[). 

We can rule out large decreases of Pzz in time between the LSP and 
LDP measurements. If the polarization were to change, it must do so in 
particular ways. If the dissociator of POLIS were to operate anomalously 
(a decrease in efficiency), then the polarization of both vector and tensor 
polarized beams would decrease. This is not observed since only the tensor 
polarized beams have lower polarization. If the medium field transitions 
units of POLIS (the transition field allowing tensor polarized beams) would 
operate anomalously, then the polarization would vary diagonally on a Pzz 
vs. Pz plot. This is not seen since, Pz of tensor polarized beams is the same 
for both the LSP and LDP, ruling out variations in the medium field transi- 
tion unit. We are left with the conclusion that the LDP underestimates the 
tensor component of beam polarization. 

The agreement between LSP and LDP can be improved if one assumes 
lower tensor analyzing powers. By reducing the tensor analyzing powers 
(A X x ~ and A zz ) of the LDP by 25%, Pzz increases as shown in Fig. 1231 
This suggest that the effective tensor analyzing powers {{A xx — A yy ) e f* and 
A e JJ) we used are wrong. Either because our model (Eq. ll7|) for effective 
tensor analyzing power is incomplete 5 , or because the analyzing powers we 
used in the model were wrong. We used the tensor analyzing powers pub- 
lished by Fletcher |24j for two reasons. Firstly because they were closer 
to the theoretical value predicted. Secondly because they were published 
for a number of energies, allowing us to estimate the energy dependence 

5 We did not include the effect of beam straggling, nor the angular efficiency profile of 
the detectors, nor the effect of material coating the front of the target. 
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Figure 23: LDP with reduced tensor analyzing powers (diamond symbol) is in 
better agreement with the LSP. 



of the analyzing powers, and use beams of any energy between 25 and 80 
keV. The tensor analyzing powers published by Becker [23] were not used 
because they were inconsistent with the theoretical predictions, and also 
because they were published at a single energy. The difference in tensor 
analyzing power between these two publications is considerable. For exam- 
ple, at Ed 30 keV and 6 cm = 90°, Fletcher j^lj has A xx — A yy and A zz 
approximately 35% and 50% greater than Becker Our effective ana- 

lyzing powers were derived from Fletcher's analyzing powers, yet our data 
would fit better with lower tensor analyzing powers, such as those published 
by Becker. We could not use Becker's tensor analyzing powers because our 
measurements was performed at different energies. 

Although Pz depends on the value of Pzz (ei m Eq. ()12[)). this depen- 
dence is very weak. Having the wrong Pzz changes Pz by a negligible 
amount. 

The vector polarization that the LDP measures is reliable because the 
vector analyzing power is energy independent, matches the theoretical pre- 
diction, and does not depend on our choice of model for effective analyzing 
powers. 

Finally, the LDP measures Pz's of vector polarized beams that are on 
average 7% lower than the LSP. This relative difference can originate from 
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unknown systematic uncertainties of the LSP or LDP, or because of polar- 
ization change in space (between the polarimeters) or in time (between the 
measurements) . 
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6 Conclusions 



We built a Low Energy Deuteron Polarimeter (LDP), based on a reaction 
that received little attention for polarimetry. We find that the LDP is well 
suited to measure asymmetries but presently lacks a proper tensor analyzing 
power calibration. 

Data from our polarization cross-check experiment cannot answer whether 
there are polarization changes during acceleration, because our IBP data are, 
at this point, preliminary. 

The relative difference, of vector polarized beams, between the LSP and 
LDP was 7% in our experiment. 

6.1 Future Improvements 

We can point out future improvements to the experiment we performed, and 
to future versions of the LDP. 

6.1.1 Polarimeter Cross-Check Experiment 

• The online analysis program of the LSP, although fast, often returns 
inaccurate polarization and uncertainties in polarization. Offline anal- 
ysis, which was done here, decreased the uncertainties (sometimes by 
a factor of 10) and changed the polarization (mostly the tensor com- 
ponent of polarization). The online analysis program could be made 
more accurate. 

• Use all the polarization states that POLIS can provide (13), not just 
pure vector (2) and tensor beams (2) but mixed vector-tensor beams 
(9), for additional systematic checks between polarimeters. 

• Our conclusions rely to a large extent on the fact that the beam polar- 
ization from the source was constant during the experiment. During 
next experiment the polarization should be monitored with the LSP 
before and after each run. 

6.1.2 LDP 

• Find better tensor analyzing powers. Either by improving the model 
that calculates effective analyzing power, or by using a better data set 
of published tensor analyzing powers. 

• Measure the energy dependence of the vector analyzing power. We 
expect it to be energy independent, but do not know to what level. 
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